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Germanium is a widely used material for electronic and optoelectronic devices and recently it has 
become an important material for spintronics and quantum computing applications. Donor spins 
in silicon have been shown to support very long coherence times {T 2 ) when the host material is 
isotopically enriched to remove any magnetic nuclei. Germanium also has non-magnetic isotopes 
so it is expected to support long T 2 S while offering some new properties. Compared to Si, Ge 
has a strong spin-orbit coupling, large electron wavefunction, high mobility, and highly anisotropic 
conduction band valleys which will all give rise to new physics. In this Letter, the first pulsed electron 
spin resonance (ESR) measurements of T 2 and the spin-lattice relaxation (Ti) times for ^®As and 
donors in natural and isotopically enriched germanium are presented. We compare samples with 
various levels of isotopic enrichment and find that spectral diffusion due to ^^Ge nuclear spins limits 
the coherence in samples with significant amounts of ^®Ge. For the most highly enriched samples, 
we hnd that T\ limits T 2 to T 2 = 2Ti. We report an anisotropy in T\ and the ensemble linewidths 
for magnetic fields oriented along different crystal axes but do not resolve any angular dependence 
to the spectral-diffusion-limited T 2 in samples with ^®Ge. 


Germanium was the original material for transistors, 
and is now being developed for the latest semiconduc¬ 
tor electronics [T]. Recently, it has become a key ma¬ 
terial for spintronics nn and quantum computing]^ 
[7] devices. Compared to silicon, donor electrons in Ge 
have higher mobility (~ 3 times) [1] , larger wavefunctions 
(6.5 nm compared to 2.5 nm), [HI IS]) stronger spin-orbit 
coupling [in], and highly anisotropic conduction band val¬ 
leys 0. Much of silicon’s success in the quantum com¬ 
puting community has hinged on the attainability of long 
coherence times (T 2 ) exceeding seconds when Si is iso¬ 
topically enriched to have no magnetic nuclei [IIHIl]- 
Germanium also has non-magnetic isotopes so it has been 
expected to support long coherence times. In this Letter, 
we report the first electron spin coherence measurements 
for donor electron spins in Ge. We find that spectral 
diffusion due to ^^Ge limits in natural Ge samples 
whereas the spin-lattice relaxation time, Ti, limits T 2 in 
isotopically enriched Ge. The longest T 2 we measured 
is T 2 = 2Ti = 1.2 ms at 350 mK in a magnetic field 
{Bq) of 0.44 T. The low-temperature Ti fits the temper¬ 
ature dependence theorized by Roth m and Hasegawa 
m which also predicts Ti oc Bq'^. This suggests that 
considerably longer coherence times are possible at lower 
fields. 

While T 2 for donors in Ge is shorter than the times 
demonstrated for Si, Ge-based qubits have some im¬ 
portant advantages. For example the larger electron 
wavefunctions relax the lithographic requirements for ex¬ 
change coupling two donors which is important for most 
donor-based quantum computing schemes |17j . This is 


advantageous considering Ge is compatible with most of 
the same nanofabrication techniques as silicon and single¬ 
donor devices are achievable m- Another useful feature 
of Ge is the large spin-orbit coupling and shallow donor 
depth which leads to a very large spin-orbit Stark shift in 
Ge (nearly 5 orders of magnitude larger than in silicon) 
[B] meaning that Ge based qubits are extremely tuneable. 
This will be important for gated quantum devices [T7]. 

Despite these features, the spin coherence of donor 
electrons in Ge has remained mostly unstudied. The 
first experiments were conducted over fifty years ago 
by Eeher, Wilson, and Gere mm: but their measure¬ 
ments were limited to continuous wave (GW) ESR spec¬ 
troscopy. They estimated Ti for ^^As and donors 
based on power saturation measurements, but experi¬ 
mental errors were large. These experiments are difficult 
because wavefunction overlap occurs for densities as low 
as 10^^ donors/cm^ such that only lightly doped sam¬ 
ples with correspondingly weak signals are useful for iso¬ 
lated donor experiments. Some limited experiments on 
Sb [HHl [S] and [55] donors in highly strained Ge 
were also reported. More recently, pulsed nuclear mag¬ 
netic resonance studies were conducted on ^^Ge nuclear 
spins [55U5H] which found that the ^^Ge nuclear spin co¬ 
herence in germanium can be > 100 ms. 

The samples discussed in this Letter include commer¬ 
cially available, natural Ge doped either 10^® As/cm^ or 
10^^ P/cm^. "^^Ge is the only naturally occurring isotope 
of Ge (7.75% abundance) with a nuclear spin and is thus 
expected to be a limiting factor in the donor spin coher¬ 
ence at low temperatures. Three isotopically enriched 


2 


samples were prepared at Lawrence Berkeley National 
Laboratory. The first is a piece of neutron transmuta¬ 
tion doped ^'^Ge described in Ref. [261 El]. This sample 
is uniformly doped with ^®As to a density of 3 x 10^® 
donors/cm^ and contains a residual 3.8% ^^Ge. The 
other two samples are 96% ™Ge crystal (0.1% ^^Ge) 
and a 99.99% ™Ge crystal (0.01% ^^Ge). They have 
concentrations of ~ 10^^ donors/cm^ and ^ 10^^ 
donors/cm^, respectively and are described in [5^1 HH]- 
The crystallographic orientation of the samples was de¬ 
termined using X-ray diffraction. The sample details are 
summarized in Table jH 

The experiments down to 1.65 K were performed in a 
pumped He cryostat (H.S. Martin), and lower temper¬ 
ature data were obtained in a ^He cryostat (Janis Re¬ 
search). All data were taken at X-band (9.65 GHz) in 
a Bruker dielectric resonator (MD5). The ESR spec¬ 
tra were measured via echo-detected field sweeps using 
a standard Hahn-echo pulse sequence (7r/2-r-7r-r- 
echo). Typical spectra are shown in Fig. [^a) for phos¬ 
phorus donors in the 0.1% ^^Ge:P sample and in Fig. Bb) 
for arsenic donors in the 3.8% ^^Ge:As sample. From 
these plots we extract a hyperfine coupling constant of 
3.55 mT for ^®As and 2.04 mT for ^^P. 

The FSR linewidth depends strongly on the sample 
orientation and the abundance of ^^Ge present in the 
sample, as noted by Wilson |8]. With Bq oriented along 
one of the (001) directions, the linewidth is narrowest and 
is limited primarily by hyperfine interactions with ^^Ge. 
At this orientation the line broadening from spin-orbit 
strain effects is suppressed by valley symmetry about the 
(001) as explained in Refs.|8l [HI [T6]L To give a sense 
of the strain-induced line broadening for Bq away from 
[001] equivalent directions, Fig0 c) shows the angular 
dependence of the linewidth for select samples rotated in 
the (110) plane relative to the [001] axis. There is also 
an isotopic dependence of the linewidth away from the 
(001) direction and we presume this is due to isotopic 
strain [30]. The strong dependence of the linewidth on 
field orientation conveniently allows for accurate in situ 
orientation of the crystals. Unless otherwise noted, all 
data presented in this manuscript assumes Bq is oriented 
along a (001) axis. 

One can predict the effect of ^^Ge on the FSR 
linewidth through the hyperfine interaction using a sec¬ 
ond moment calculation m, which gives AB oc 
where AB is the linewidth, and / is the percent abun¬ 
dance of ^^Ge. The measured FSR linewidths for samples 
of various isotopic enrichment with i?oll(001) is shown in 
Fig. T[d). The point at / = 0.8% was taken from Wil¬ 
son ^ . The solid curve in Fig. B'l gives the expected 
dependence for broadening of the line due to ^^Ge 
hyperfine interactions for ^^As. The solid curve fits the 
data well, implying that ^^Ge is indeed the dominant 
mechanism for line broadening in this orientation. The 
linewidth can be interpreted as an ensemble dephasing 




FIG. 1. (a) Echo-detected field sweep spectra for (a) 0.1% 

^®Ge:P and (b) 3.8% '^^Ge:As with Bo || (001). (c) Plot of 
ESR linewidths as a function of field orientation for "“*Ge;As 
(blue), 3.8% ^®Ge:As (red) and 0.01% ^^Ge:P(black). The 
solid lines serve only as guides to the eye. (d) Linewidth for 
Bo II (001) as a function of ^^Ge isotopic abundance. The 
Ge:As data appear as black triangles whereas the Ge:P data 
appear as red circles. The solid line shows the expected 
dependence for broadening due to ^^Ge hyperfine interactions. 
The ESR linewidth at 0.8% is taken from Ref. [Hj. Data were 
taken at 1.8 K and 9.65 GHz. 


time, T 2 , which is also shown in Table [T| 

Ti was measured using an inversion-recovery pulse se¬ 
quence (7r-t-7r/2-r-7r-r - echo). The values of 
Ti are plotted in Fig. for ^^P(a) and ^®As(b) donors. 
The same two mechanisms limit Ti for all of the sam¬ 
ples. At higher temperatures, Ti is limited by a highly 
temperature (T) dependent process. The theory of Roth 
and Hasegawa [HI [16] predicted a Raman process 
to dominate at these temperatures but this dependence 
does not fit our data well. An Orbach process does fit 
the data as shown in Fig. [^ The Orbach process is of 
the form Ti oc a x exp{ ), where a is a prefactor that 
can be calculated using Ref. [32], F^.o. is the valley-orbit 
splitting, and k is the Boltzmann constant. The valley- 
orbit splittings extracted from the Ti fits in Fig. [^ agree 
well with the values measured by Ramdas (2.8 meV for 
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TABLE I. Sample Details 


Sample Name 

[^°Ge] 

[""Ge] 

[""Ge] 

[""Ge] 

[""Ge] 

Doping (cm ®) 

[001] Linewidth (mT) 

TI (ns) 

"“*Ge:As* 

20.57% 27.45% 

7.75% 36.50% 

7.73% 

1 X 10^® As 

1.2 

11 

"“‘Ge:P* 

20.57% 27.45% 

7.75% 36.50% 

7.73% 

~ 10^2 P 

1.1 

13 

3.8%^3Ge:As 

0.1% 

0.9% 

3.8% 

92.6% 

2.6% 

3 X 10^® As 

0.8 

17 

0.1% ''^Ge:P 

96.3% 

2.1% 

0.1% 

1.2% 

0.3% 

~ 10^^ P 

0.069 

211 

0.01% ’’3Ge:P 

99.99% 

- 

0.01% 

- 

- 

~ 10 “ P 

0.051 

284 

Nuclear Spin 

0 

0 

9/2 

0 

0 

- 

- 

- 


‘Percent abundances for the natural germanium samples were taken from Ref. 




FIG. 2. Temperature dependence of T\ (triangle) and T 2 
(circle) for natural (open symbols) and isotopically-enriched 
(solid symbols) Ge with Bq || (001). The solid lines are fits for 
(a) phosphorus donors (0.1% ^®Ge:P) and (b) arsenic donors 
(3.8% ^®Ge:As), assuming two relaxation processes: a single¬ 
phonon {T~^) process and an Orbach (a x exp{ )) process. 
For the T 2 fits, both T\ and an additional (temperature in¬ 
dependent) spectral diffusion mechanism due to ^^Ge were 
taken into account. Note that for the 0.1% ^®Ge:P sample, 
T 2 = 2Ti down to the lowest temperatures. 


and 4.2 meV for [H]). Likewise, the values of a 
extracted from our fits agree with the values calculated 
using Castner’s theory [32] to within a factor of 2. 

At lower temperatures, a single-phonon process with 
a T~^ dependence appears to dominate. This relaxation 
process is a result of the multivalley structure of germa¬ 
nium. In the unperturbed ground state, there are four de¬ 
generate valleys located along the (111) equivalent crys¬ 
tallographic axes. Each valley has an axially symmetric 
g-tensor, gi but the effective g-tensor, geff, is given as 
a weighted average over all four valleys. In the electron 
ground state, each valley has equal amplitude, and, by 
symmetry, geff is isotropic |15j . When strain is applied, 
valley energies shift relative to each other, leading to val¬ 


ley repopulation and a change in *gef /. The strain from 
phonons near the Larmor frequency modulates geff^ ef¬ 
fectively mixing the spin up and down states. This gives 
a Ti as calculated by Roth m and Hasegawa m which 
agrees well with our experimental data. The calculated 
estimates for Ti at 350 mK are within 10% for Ge:As 
and 30% for Ge:P. The theory predicts that Ti due to 
this single-phonon process should scale with the square of 
the *gi anisotropy. The valley anisotropy of Ge was mea¬ 
sured to be 3 orders of magnitude larger than in Si |H], 
implying that the single-phonon process should be 6 or¬ 
ders of magnitude stronger in germanium. This accounts 
for the short Ti times observed for donors in germanium 
as compared with silicon. 

An interesting property of the single-phonon spin- 
lattice relaxation mechanism is an anisotropy in Ti pre¬ 
dicted by the Roth-Hasegawa theory[T5l [16]. The 3.8% 
'^^Ge:As crystal was rotated in the (110) plane at 1.8 
K, and the resulting Ti is plotted in Fig|^ The theory 
predicts that, for rotation in this plane, the spin-lattice 
relaxation is given by: 

^ = aB*T{cos'^{9) + ^sin'^{e)) (1) 

where a is a scaling factor which can be calculated follow¬ 
ing Hasegawa|16]. and 9 is the field orientation relative 
to (001). Hasegawa calculated a = 7.2 x 10^K~^s~^T~^ 
for arsenic in Ge, but a fit to the data reveals a = 
4.1 X . We observe that for Bq oriented 

along a (111) axis, Ti becomes 3 times longer than along 

( 001 ). 

We note that Ti for donors in highly enriched samples 
is shorter than it is for donors in the natural material as 
seen in Fig. [^b). This effect is still under investigation, 
but one possible mechanism is the presence of isotopic 
strain in the natural germanium [30]. Wilson [8] demon¬ 
strated the use of large strains to partially lift the valley 
degeneracy, thus disrupting the single-phonon relaxation 
mechanism. Modelling the effects of strain can be com¬ 
plex, as strain not only modulates a, but can also modify 
the form of Eq.([^. Nevertheless, controlled strain may 
be beneficial for future quantum devices based on germa¬ 
nium. 

We also measured the electron spin coherence time. 
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FIG. 3. Angular dependence of Ti for the 3.8% ^®Ge:As 
sample rotated in the (110) plane at 1.8 K. The curve is a fit 
using Eq.Q assuming a = 4.1 x 10'* K~^. 

T 2 , for each of the samples using the standard Hahn-echo 
pulse sequence. The decay curves at 1.8 K for Bq || (001) 
are shown in Fig.|^a) for Ge:P and in Fig.|^b) for Ge:As. 
These decays are fit to an exponential decay of the form 
^ where A scales the amplitude, r is the delay 
between the 7r/2 and tt pulses in the Hahn echo sequence, 
and n is a fitting parameter that depends on the decoher¬ 
ence mechanism. The 0.1% ^^Ge:P sample decays with 
ri=l over the measured temperature range. For this sam¬ 
ple it was found that T 2 = 2Ti (representing the absolute 
Ti limit [34] ) down to 350 mK temperatures, meaning 
that decoherence due to ^^Ge is negligibly small with 
this level of isotopic enrichment at these temperatures. 
For samples with / > 3.8%, we find that n varies from 
1 at high temperatures to 2.1 at low temperatures. This 
is a characteristic of ^^Ge spectral diffusion limiting the 
coherence. At 1.8 K, the "“‘Ge:As, "“*Ge:P, and 3.8% 
^^Ge:As samples decay with this form. 

The temperature dependence of T 2 is also plotted in 
Fig. and fit to I/T 2 = l/(2Ti) -|- I/Tsd, where Tsd is 
the (temperature independent) spectral-diffusion-limited 
coherence time. For the natural germanium samples, 
Tsd limits the coherence to 57 /rs whereas the 3.8% 
^^Ge:As sample is limited to 113 /iS. From similar work 
in silicon [3SJ|3n]j one might expect an orientation depen¬ 
dence to Tsd- We measured the orientation dependence 
of Tsd for the 3.8% ^^Ge:As sample at 1.8 K and fit the 
decays with a curve of the form 

to separate the Ti component from Tsd 137]. No angular 
dependence of Tsd could be resolved. 

While coherence times of over one millisecond for iso- 
topically enriched material open the possibility of using 
donor electrons in Ge for quantum computing devices, 
these coherence times are much shorter than those for 
donors in isotopically enriched silicon (seconds) [TU [T3| . 
To extend the Ge donor coherence, one must either over¬ 
come the Ti limit or use nuclear spins which may sup¬ 
port longer coherence times. There are several promis¬ 
ing techniques to extend the Ti limit. One approach is 
to take advantage of the Ti anisotropy, which will al- 



FIG. 4. Two-pulse Hahn echo decay curves for natural 
(blue) and isotopically enriched (black) germanium doped 
with phosphorus(a) and arsenic(b) donors. Data were taken 
at 1.8 K and 9.65 GHz. The solid curves are fits to the data 
using exp[-{^)’^]. 


low for up to a factor of 3 increase in T) when devices 
are oriented with Bq || (111), but this Ti enhancement 
comes at the expense of a shorter ensemble T 2 . A simple 
alternative is to operate devices at lower temperatures, 
since Ti oc T“*. Perhaps the most effective technique is 
to operate devices at lower frequencies since theory pre¬ 
dicts Ti oc Bq'^. More complicated strategies are also 
available. In particular, one can apply a large strain, as 
demonstrated by Wilson [5] which shifts the valley energy 
levels, thus suppressing valley repopulation and the asso¬ 
ciated relaxation mechanisms. Another recent proposal 
suggests patterning Ge in a periodic structure to open a 
phononic bandgap at the Larmor frequency [38]. Such a 
structure would suppress the single phonon process. 

In summary, we have measured the ESR linewidths, 
coherence times, and spin-lattice relaxation times for 
donors in natural and isotopically enriched germanium 
at X-band microwave frequencies. We find that the 
linewidths are primarily broadened by hyperfine inter¬ 
actions with ^^Ge spins when Bq is oriented along the 
[001] axis and by strain in other orientations. We find 
that donor electron spin coherence is limited by spec¬ 
tral diffusion due to hyperfine interactions with ^^Ge 
nuclei for the "“*Ge {Tsd = 57 p,s) and 3.8% "^^GeiAs 
{Tsd = 113 pis) samples, thus Tsd scales approximately 
as 1// which is similar to silicon|3B]. For the more 
highly enriched 0.1% ’^^Ge:P sample, T 2 was limited to 
2ri down to 350 mK, the lowest temperature we have 
measured (T 2 = 1.2 ms for Bq || (001)). We observe a 
large anisotropy in Ti, which is explained by the theory 
of Roth and HaseeawafTOl ITS] , with the longest Ti occur- 
ing for Bq || (111). It is predicted that at lower magnetic 
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fields Ti and thus T 2 should become substantially longer. 
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